Odd integer quantum Hall effect in graphene 
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A possible realization of Hall conductivity, quantized at odd integer factors of e^ /h for graphene's 
honeycomb lattice is proposed. I argue that, in the presence of uniform real and pseudo-magnetic 
fields, the valley degeneracy from the higher Landau levels can be removed. A pseudo-magnetic 
field may arise from bulging or stretching of the graphene flake. This may lead to observation 
of plateaus in the Hall conductivity at quantized values fe^/h, with / = ±3, ±5 etc, which have 
not been observed in measurement of Hall conductivity. However, in a collection of noninteracting 
Dirac fermions living in the honeycomb lattice subject to real and pseudo field, the zeroth Landau 
level still enjoys the valley and the spin degeneracy. Upon including the Zeeman coupling, the spin 
degeneracy is removed from all the Landau levels. The effects of short ranged electron-electron 
interactions are also considered, particularly, the onsite Hubbard repulsion (U) and the nearest- 
neighbor Coulomb repulsion (V). Within the framework of the extended Hubbard model with only 
those two components of finite ranged Coulomb repulsion, it is shown that infinitesimally weak 
interactions can place the system in a gapped insulating phase by developing a ferrimegnatic order, 
if i7 >> y. Therefore, one may expect to see the plateaus in the Hall conductivity at all the integer 
values, / = 0, ±1, ±2, ±3, • • •. Scaling behavior of interaction induced gap at / = 1 in presence of 
finite pseudo fiux is also addressed. Qualitative discussion on finite size effects and behavior of the 
interaction induced gap when the restriction on uniformity of the fields are relaxed, is presented as 
well. Possible experimental set up that can test relevance of our theory has been proposed. 

PACS numbers: 71.10.Pm, 71.70.Di, 73.43.Nq 
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I. INTRODUCTION 

Carbon atoms bonded in a two-dimensional honey- 
comb lattice, named as graphene recently engaged 
both theoretical and experimental attention due to a va- 
riety of unconventional properties. The lack of inversion 
symmetry in graphene lattice structure leaves the Bril- 
louin zone with special points at its corners, allowing 
one to linearize the spectrum around those points, which 
leads to a conical dispersion. At filling one-half, when the 
chemical potential coincides with the apex of the cone, 
excitations can be described in terms of massless, chiral 
Dirac fermions, where the Fermi velocity (i;^ = C/300) 
plays the role of the velocity of light (C) [4]. 

Many peculiarities in graphene's electronic properties 
arise from the Dirac nature of the quasi- particles. One 
of the examples is the anomalous quantum Hall effect. 
Placed in a magnetic field graphene exhibits integer 
quantum Hall effect, where the Hall states reside at 
filling factors u — ±(4n -I- 2) with n — 0,1, 2,-- - at 
low magnetic fields 10 T ) [3]. This phenomena 
can be understood from the relativistic nature of the 
non-interacting quasi-particles. Upon exposing the 
system to stronger magnetic fields (> 20 T), additional 
Hall states appear at filling factors v = 0,±1,±4 
However, the plateaus in the Hall conductivity at other 
odd integer fillings ( for example, v = ±3 or ±5 etc. 
) are absent even at the highest laboratory magnetic 
field 45 T ) 5]. These observations confirm that the 
four fold degeneracy of the zeroth Landau level (LL) is 
completely lifted. Furthermore, the sub-linear variation 
of the activation gap with magnetic field of the ly — 1 
Hall state points its origin toward electron-electron 



interactions. Whereas the linear dependence of the 
activation gap with total magnetic field for v — ±4 Hall 
states suggests that, possibly a finite Zeeman coupling 
removes only the spin degeneracy from the higher LLs. 
Within the framework of the extended Hubbard model 
with onsite and nearest-neighbor repulsion, one can 
show that the valley degeneracy only from the zeroth 
LL can be lifted, due to a spontaneous breaking of 
the chiral symmetry. In presence of a magnetic field 
such metal-insulator transition can take place even at 
sufficiently weak interaction. However, the higher LLs 
gain a finite shift in energy only. Therefore, one can 
explain the absence of Hall states at odd integer fillings ( 
e.g., 1/ = ±3, ±5 ) from the protected valley degeneracy 
of the higher LLs [1,0]. 



In addition to the two-dimensional crystal struc- 
ture, ripples are also present on the graphene sheet 
and introduce a fictitious gauge potential in the long 
wavelength limit 8]. Such corrugation of the graphene 
sheet is partially due to the strain induced by the 5^02 
substrate Q. In a normal graphene flake these defects 
are randomly distributed and therefore give a net zero 
flux of the fictitious field. Perhaps, bulging or stretching 
of the graphene sheet can give rise to a finite pseudo fiux 
. Recently, its been argued that specific distortions of 
a fiake can give rise to a uniform pseudo-magnetic field 
up to 10 T [11]. Nevertheless, a pseudo-magnetic field of 
strength ~ 350 T, is produced by depositing a graphene 
layer on platinum substrate, followed by cooling of the 
system [12]. In contrary to the real magnetic field, the 
pseudo vector potential couples with opposite sign at 
two Dirac points and hence preserves the time reversal 
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symmetry (TRS). Therefore, in the presence of uniform 
real and pseudo field, one expects the LLs to loose the 
valley degeneracy. This may also allow the formation 
of plateaus in Hall conductivity at odd integer values of 
/h which were previously absent. 

Our discussion is concerned with quantization of 
the Hall conductivity in graphene in the presence of 
both real and pseudo magnetic fields. The rest of the 
discussion is organized as follows. In Sec. H, we compute 
the spectrum of the free Dirac Hamiltonian in the pres- 
ence of both fields. The short ranged electron-electron 
interaction is discussed in Sec. HI. Section IV is devoted 
to the scaling behavior of interaction induced gap. We 
present the concluding remarks and discuss some related 
issues in Sec. V. 



II. FREE ELECTRON SPECTRUM 



after conveniently rotating the frame of reference to 
qx — q - KjK and qy = [K x q) x KjK"^ . We adopt 
natural units, h = e — vp — \, where vp — ta\/i/2 is the 
Fermi velocity. The cut-off in the momentum integral 
( A ^ 1/a ), represents the interval of the energy over 
which the tight binding density of states is approximately 
linear. The Einstein summation convention is assumed, 
but only over repeated space-time indices. The mutually 
anti-commuting four component Hermitian gamma 
matrices belong to the "graphene representation", 
lo = I2® era, 71 = CTa 18) CT2 and 72 = I2® cri. We define 
the remaining two gamma matrices as 73 = cri ® (T2 and 
75 = a2® (J2- The low energy Hamiltonian Hq preserves 
the emergent "chiral" C/c(4) symmetry generated by 
{jSbp® {^4,73,75735} where, 735 = 17375 = a^, ® I2 
[l5|, [l6|. The two component Pauli matrices {to,t}, 
operate on spin indices. One can study the response 
of magnetic field by defining Di = —idi — Ai, with 
magnetic field B = e^ijdiAj set to be perpendicular to 
the graphene plane. 



The tight binding Hamiltonian for spin- 1/2 electrons 
on the honeycomb lattice in the presence of only nearest- 
neighbor hopping is defined as 



(1) 



Here u^(^) denotes the fermionic creation operator on 
the sites of one of the triangular sublattices generated 
by the linear combination of the basis vectors ai = 
(•\/3, — l)a and 02 = (0, l)a. a is the lattice spacing ( 
^ 2.5 A ) and t is the nearest-neighbor hopping am- 
plitude ~ 2.5 eV (13i]. Analogously, v{A + bi) is the 
electron annihilation operator on the other sublattice, 
then located at i? = A + b, with the vector b be- 
ing either 61 = (l/\/3, l)a/2, 62 = (1/a/3, -l)a/2 or 
&3 = (— 1/1/3, 0)a. Within the framework of the tight 
binding model, energy spectrum E = ±t\ exp {ik ■ bi)\ 
is doubly degenerate and becomes linear and isotropic 
in the vicinity of the corners of the first Brillouin zone. 
Among six such points, only two are inequivalent and 
suitably chosen to be at ±K, with K = {1, l/-\/3)27r/a\/3 
14{ . Hence, retaining the Fourier modes only near these 
two points, one can write down the effective Hamilto- 
nian corresponding to the tight binding model in the low 
energy limit as, 



(2) 



with 

^l{x, t) 



A 



dq 



;e-'^-%l{K + q)A{K + q), 



(27ra)2 

ul{-K + q),vl{-K + q)], 



In addition to the Dirac quasi-particles, a crucial in- 
gredient of the system is ripples. The effect of the rip- 
ples can be captured in terms of fictitious gauge fields 
in the long wavelength limit. In graphene fiakes, ripples 
are randomly distributed and therefore give a net zero 
fiux of pseudo-magnetic field. Whereas, a bulging of the 
graphene sheet may introduce a finite flux of the pseudo- 
magnetic field. The Dirac Hamiltonian in presence of 
both the real and pseudo-magnetic fields reads as 



Ho[A, a]^TQ® ijoliiPi - Ai - aj'^735). 



(4) 



where, ap is the member of a general non-Abelian SU{2) 
gauge field. 



a-i = a,73 0^75 + 735. 



(5) 



(3) 



A smooth enough deformation in the graphene sheet, do 
not mix two inequivalent valleys at K and —K. Assum- 
ing that the bump in the graphene flake varies slowly on 
the lattice scale we kept only one component of a general 
SU{2) gauge potential [iTj. One way to introduce 
both gauge potentials in experiment is the following. 
First deposit graphene over a metallic substrate, e.g., 
platinum [l^ , at a relatively high temperature and then 
cool the system. Due to a mismatch in compressibility 
of the substrate and graphene, the former one produces 
strain on the graphene sheet. This way one might expose 
the system first to a finite pseudo-magnetic field. Once 
the fictitious field is introduced, the system can then be 
placed in a real magnetic field. 

It is informative to cast the Hamiltonian _ffo[^, a] in 
the following block-diagonal form, Ho[A, a] = H^[A, a] © 
H^[A, a], where 

H± = ±l2<S)(Ji{-idi-AiTaf)-l2<»cr2{-id2-A2Taf)- 

(6) 
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H+ and H- represent the Hamiltonian near K and —K 
point, respectively. From Eq. [B] one can register that the 
net magnetic field near one Dirac point [K) is enhanced 
to the value B^^^^i = B + b, with b = ezijdiO^^ . Near 

the other Dirac point at —K the net magnetic field is 
attenuated to B^^^^i — B—b. Unless mentioned otherwise 
we assume S > 6 for the rest of our discussion. Both the 
fields are assumed to be uniform as well. One can see 
that both H± are unitarily equivalent to a generic Dirac 
Hamiltonian H£,[A^ a] in two dimensions in the presence 
of gauge fields, where 

Hd[A, a] = i7o7*(Pi -Ai-af). (7) 

Specifically, = ulHD[A,a]Ui, with Ui = h ® icF2 
and = uIHd[A, -a]U2, with U2 = i(72 ® h [3- 

In the absence of the pseudo-magnetic field (6 = 0) 
the Hamiltonian i7o[^, 0], exhibits a series of LLs at 
well separated energies ±\'2nB, n = 0,1,2,---, with 
degeneracies oi B/tt per unit area. At half- filling, all the 
negative energy LLs are completely filled, whereas the 
LLs at positive energies are totally empty, with only half 
of the zero energy states occupied. As a consequence, 
the Hall conductivity shows plateaus at integer fillings 
1/ = ±(4n -f 2), n = 0,1, 2,--- Measurement of 

the Hall conductance at relatively low magnetic fields 
{B ~ 10 T) confirmed such quantization [3, 4]. The 
additional four-fold degeneracy of the LLs arises from 
the spin and the valley degrees of freedom. The same 
story follows when one switches off the real magnetic 
field, and only a pseudo-magnetic field penetrates the 
system, with the difference that now the LLs appear at 
energies ±^/2nb. A recent experiment confirmed such 
quantization in the presence of a pseudo-magnetic field 
of strength ^ 350 T J^]. However, the presence of finite 
real and pseudo-magnetic field removes the valley degen- 
eracy fro m ea ch LL. Consequently, each of the LLs at 
energy upon imposing a finite pseudo-magnetic 

field, gives rise to two new LLs at energies y/2n{B ± b), 
with degeneracies D± — {B zL &)/27r per unit area, 
respectively, but only for n 7^ 0. It is worth noticing 
that the Dirac points are decoupled from each other 
even in presence of a finite pseudo field. Hence, there 
are two sets of zero energy states, one is localized near 
the K point, and the other one near —K (or K') with 
degeneracies $7(5 ± 6)/27r, respectively. Here 17 being 
the area of the sample. Thus, the zeroth LL does not 
split even when the system experiences finite real and 
pseudo fields [20I [2l| . However, in presence of uniform 
real and pseudo- magnetic field the higher LLs split, 
thereby pushing the states near K point up in energy, 
and those in the vicinity of K' point down in energy. 
Hence, with some particular strength of B and 6, states 
from two successive LLs, localized near different Dirac 
points can be degenerate. This situation can easily be 
bypassed by tuning the ratio B/b close to an even integer. 



Upon including a finite Zeeman coupling, the spin 
degeneracy can be lifted from all the LLs including 
the zeroth one. The Zeeman splitting scales as Az (in 
Kelvin) ~ B, where B is measured in Tesla, thus much 
smaller than the LL energies. For our purposes we tune 
the magnetic fields so that, 5 < B/b < 10. Within 
this range of parameters, the energy difference between 
the Hall states with same LL index (n) but localized 
near two different valleys is ~ 200 — 400 K. Energies 
corresponding to different LLs for one particular set of 
realistic values of the magnetic fields are shown in Table 
I. Without considering the many body effects arising 
from the electron-electron interactions, one can expect 
the Hall conductivity to exhibit plateaus at integer 
values / = 0,±2,±3,±4,±5,--- of e'^/h. It is worth 
noticing that the LLs associated with different Dirac 
points do not enjoy the luxury of equal degeneracy. 
Therefore, plateaus of Hall conductivity at integer values 
of e'^/h, are not associated with the integer fillings. 
For further illustration of the computation of the Hall 
conductivity the reader may refer to the Appendix. 



TABLE I: Energies of LLs with B = 32 T and fe = 4 T 
n 1 2 3 



Eo 2.4 X 10^ K 3.4 x 10^ K 4.1 x 10^ K 
E+ 2.5 x 10^ K 3.5 x 10^ K 4.4 x 10^ K 
E- 2.2 X 10^ K 3.2 x 10^ K 3.8 x 10^ K 



Energies of the LLs (measured from the charge neutral point) for 
a particular choice of the magnetic fields. Here n stands for the 
LL index in absence of pseudo field. Eq stands for the energy of 

the LL without any pseudo field. -E+(_) denotes the energy of the 
LL lives in the vicinity of the +K(~K) Dirac point. 



III. ELECTRON-ELECTRON INTERACTIONS 

Diagonalizing the Dirac Hamiltonian in the presence 
of uniform real and pseudo-magnetic fields we found that 
the valley degeneracy is lifted only from higher LLs. In 
the absence of interactions the zeroth LL still enjoys the 
valley degeneracy. In this section we will consider the 
effect of the short-range electron-electron interactions in 
the spectrum. The interacting Hamiltonian in presence 
of only on-site (U) and nearest-neighbor (V) repulsion is 
defined as, 

-ffc/ = + Y n„{A)n^,{A + bl). 

(8) 
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For graphene, /7 « 5 - 12 eV and U/V « 2-3 [12]. 
We consider the system to be at filling one half. In the 
absence of magnetic fields, a sufficiently large on-site in- 
teraction can take the system into an insulating ground 
state with the magnetization altering its sign at each 
site from its neighbor. A staggered pattern in average 
electron density can be realized at large enough nearest- 
neighbor Coulomb repulsion. The transitions out of the 
semi-metal to the Mott insulators are believed to be con- 
tinuous and belong to the Gross-Neveu universality class 
[isl [iGj . Therefore, we assume a uniform background of 
either the staggered density [C — {n{A) — n{A + h))] or 
staggered magnetization [N = {m.{A) — m{A + h))]. Here 
n{A) = u\{A)u„{A) is the average electron density and 
m{A) = ul{A)T'iUa{A) corresponds to average magneti- 
zation on sublattice A. Similar quantities are analogously 
defined on B sublattice in terms of fermionic operators 
Va{B) and vl.{B). After the usual Hartree-Fock decom- 
position, one can write down an effective single-particle 
Hamiltonian in that background as 

Hhf = To<g) Ho[A,a] + maig) ^0, (9) 

where a = tq corresponds to m = C [ charge density 
wave (CDW)] and a = T3 to m = [ anti-ferromagnet 
(AF)]. Here we omitted the Zeeman term for conve- 
nience. The effect of Zeeman splitting will be discussed 
later. 

The spectrum of the Hamiltonian Hjjp is as follow: 
For n 7^ the eigenvalues are at ±y/2n{B ± 6) -I- m^, 
for each spin projection, with degeneracies per unit 
area D± = {B ± 6)/27r, respectively. Besides these, for 
n = 0, the eigenvalues of Hhf for each spin components 
are ±|m| with degeneracies D° — {B ± b) /2n per unit 
area, respectively. Therefore filling up only the states at 
negative energy while leaving those at positive energy 
empty one immediately develops a gap even at small 
enough interaction. Generating a gap at infinitesimal 
interactions in the presence of a magnetic field is termed 
as "magnetic catalysis" [2^. This has been proposed 
as a mechanism behind formation of the Hall states 
in graphene at filling factor v = and 1 in presence 
of an ordinary magnetic field @, 0. In the zeroth 
LL, states associated with K and —K are localized 
on complimentary sublattices A and B respectively. 
Hence, for each spin projection, there are {B + b)/2n 
states per unit area on the A sublattice, whereas the 
other sublattice hosts only {B — 6)/27r states. In the 
absence of pseudo-field, states on each of the sublattices 
enjoy equal degeneracy. Therefore, even infinitesimally 
strong interactions can develop a gap, at filling one-half 
by spontaneously breaking the chiral symmetry of 
the Dirac quasi-particles. However, the nature of the 
insulating ground state depends on the strength of the 
interactions at the lattice scale. For example, if [/ ^ T^, 
AF ordering lowers the the energy of the ground state. 
A CDW order can develop within the zeroth LL if on 



the other hand, the nearest-neighbor component of the 
finite range Coulomb repulsion is the dominant one. 
However, due to imbalance in the number of states 
living on two sublattices in the zeroth LL, only the 
ferrimagnetic order can be developed while keeping 
the system at charge neutrality. The electron-electron 
interactions remove either the valley degeneracy or 
the sublattice degeneracy from the zeroth LL. These 
two are equivalent however only within the zeroth LL. 
Alternatively, a large enough Zeeman coupling will 
orient the magnetization at all the sites in the direction 
of the magnetic field and hence lower the energy of the 
filled Dirac-Fcrmi sea. Such a ground state with finite 
magnetization can also be stabilized by increasing the 
component of the magnetic field parallel to the graphene 
plane. The exact nature of the ground state at = in 
graphene in presence of only real field is still a subject 
of debate [24| - |27| . The activation gap of longitudinal 
resistivity for 1/ = 1 was found to be independent of the 
component of the magnetic field parallel to the graphene 
plane, and varies sublinearly with the its perpendicular 
component. This observation strongly suggests that it 
is originating from the electron-electron interactions [28| . 

Therefore upon including the Zeeman splitting and 
the finite ranged interactions the four fold degeneracy 
from all the LLs is completely removed. Within the 
zeroth LL, a gap develops in by lifting the valley or 
sublattice degeneracy. Under this circumstance, one 
can expect the Hall conductivity to exhibit quantized 
jjlateaus at cTxy = fe^/h, with / = 0, ±1, ±2, ±3, • • • 



IV. SCALING OF INTERACTION INDUCED 
GAP 

Let us now focus on the scaling behavior of the in- 
teraction induced gap. Validity of our theory relies on 
the assumption that the magnetic fields are low in com- 
parison to the characteristics lattice magnetic field, or 
equivalently the magnetic length is much larger than the 
lattice spacing. This condition is easily satisfied even 
at the highest laboratory magnetic field B ~ 45T and 
< h/B < 0.2. In a current experimental situation with 
b ^ 350 T [12], the magnetic length (w 35A) is more 
than an order magnitude larger than the lattice spacing 
(a ~ 2. 5 A). Hereafter we assume that the spin degener- 
acy is completely lifted by the Zeeman splitting and the 
chemical potential lies close to the Zeeman shifted 'Dirac' 
point, thus / = 1. After integrating out the fast Fourier 
modes within the momentum shell 1/a < k < 1/Ib, one 
can write down the effective low energy Lagrangian cor- 
responding to Hu in Eq. [5] as 
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92 



(10) 



where = ^'^TOi = —idr, and r is the imag- 
inary time, /i = 0,1,2 runs over the space-time in- 
dices and summation over repeated indices is assumed. 



Here, Ib ~ is the magnetic length and we kept only 
the least irrelevant couplings gi = {iV — U)o? /% and 
(72 = Ua^ /8. Assuming a uniform background of ei- 
ther staggered electron density or ferrimagnetic order, 
the ground state energy for the N number of four com- 
ponent fermions in the magnetic fields is given by 



J 



E{m) - E{0) _ m2 
N ~ Jg 



E 



B 



47r3/2 



ds 
372 



X(sA2)(coth(s(S + ab)) - 1)], 



(11) 



r 



as iV — > cx). Setting 6 = 0, one recovers the ground state 
energy in the presence of only real magnetic field Q- 
Here K{x) is the cut-off function introduced to sum over 
the n ^ LLs. This function satisfies K{x 00) = 1 
and K(x — )• 0) = 0, but otherwise arbitrary. E{m) can 
be understood as the Hartree-Fock variational ground 
state energy of the electrons either in a CDW or ferri- 
magnetic background. Due to the different degeneracies 
of the states localized near K and K', the energy of the 
ground state is maximally lowered by pushing down all 
the states on A sublattice below the chemical potential 
while leaving those on B sublattice empty. In the absence 
of pseudo flux, A and B sublattice hosts equal number 
of states. System then spontaneously chooses the ground 
state by breaking the Ising like symmetry of either sub- 
lattice or the valley degrees of freedom. The first term in 
the parentheses in Eq. |11| counts the contribution from 



the zeroth LL whereas, the second one includes the con- 
tributions from higher LLs. Minimizing i?(rn), one can 
cast the gap equation in the following form 



2 m 



[12) 



where X — {B + b)/m? , 5 — (gA) ^ — (ffc^) ^ measures 
the deviation from the critical interaction (gc) and q = 
{B + b)/{B-b), yielding B/b= {q + l)/{q-l). The non- 
universal value of the critical interaction is determined 

by 



1 

9c 



A 



ds 



m 
,3/2 ■ 



and the function /(A, q) is defined as 



(13) 



J 



ds fsXBn 
s3/2 XB + b 



Xse- 



^2Xs 



- 1 



Xse- 



^2Xs/q _ I 



(14) 



In the limit of the low magnetic fields B,b Bq, where 
Bq ^ is the magnetic field corresponding to the lat- 
tice scale, one can substitute K{x) by unity. Upon set- 
ting (? = 1, one gets the gap equation in presence of the 
ordinary magnetic field, obtained previously [soj . How- 
ever, for 1 < g < 1.5, or equivalently > b/B > 0.2, the 
solution of the gap equation lies in the range 17.913 > 
A = Xq > 14.0, when 5 = 0. Realizing that the solu- 
tion of the gap equation for S > or g < gc will exist 
at A > Aq, one can expand /(A, q) for large A, at least 



when q is not far from unity. Keeping the terms up to 
the second order one can write 



fix,q) = « (^1 + -L^ ^^^^^ + v^)+oix-'/'), 

(15) 

with u = 1.03258 and v = 0.461808. Hence one can cast 
the self consistent equation of the gap in terms of an 
algebraic equation 



J 



_2q_ 

g + 1 



m + V 



r 



-^V5 = 0. 

1 + 9 



(16) 



At criticality, i.e., S — 0, the interaction- induced mass wise it is a universal number. Particularly for q = 1, i.e., 
varies as m = ^/2B/C, where C now depends on q, other- 
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FIG. 1: (Color online) Gap (m) at / = 1 as function of 
B, with q = 1.285 [red (thick)], i.e., B/b = 8 and g = 1 
(black), i.e., h = Q. Here m and B are measured in units 
of ufA and Bo ~ A^, respectively. The top curves corre- 
spond to the critical point 5 = and the remaining ones 
to 5 = 0.03,0.07,0.14,0.31,0.7 (top to bottom), where 5 = 
idc - g)/gcgf^- 



6 = 0, this ratio was found to be C = 5.985 + 0{1/N) 
[30|. Recently the same universal number has been 
computed numerically on a finite honeycomb lattice, 
where we found this number to be extremely close to 
the one computed in continuum limit [sij. From the 
solution of the gap equation at S — 0, one finds that 
the size of the gap at / = 1 increases upon introducing 
a finite pseudo flux. This issue can be resolved by 
considering the zeroth LL only. In the presence of real 
and pseudo fields. Hall conductivity develops a plateau 
at / = 1 by filling il{B + 6)/27r states from the charge 
neutrality, where is the area of the sample. Hence 
the zeroth LL contribution to the ground state energy is 
proportional to ^lm{B-\-h)/2'K, where m is the activation 
gap for / = 1. Therefore, one finds that the zeroth 
LL contribution is enhanced in the presence of pseudo 
flux. In the gap equation this contribution is dominant 
whereas, the higher LL contributions are 0{m^ / B±) 
and hence sub-dominant. After some tedious, otherwise 
straight-forward calculations it can be shown that the 
higher LL contribution also increases the size of the gap 
at / = 1. We computed the gap as a function of B for 
q = 1 and 1.285 for various S and a family of such curves 
is shown in Fig. 1. The first case, q = 1 corresponds to 
& = and g = 1.285 implies B/b « 8. 

The gap at / = 1 is generated by breaking either 
sublattice or the valley degeneracy within the zeroth 
LL. The semimetal-insulator transition belongs to the 
Gross-Neveu universality class. Therefore, near transi- 
tion the dynamical critical exponent z is equal to unity 
and the correlation length exponent 1/^1 + 0{1/N) 
[l6j . The 1/A^ corrections are found to be small in 
comparison to the leading order contribution. This 
is the reason behind their omission here (32| . One 
can also confirm the critical exponents from the linear 



variation of the gap with magnetic field B, for small 
interactions {S <^ 0) As one enters the regime of 

stronger interaction, a sub-linear dependence of the gap 
on magnetic field emerges. At criticality (5 = 0, the gap 
shows a perfect square-root dependence on magnetic 
field [lOl- For q — 1.285, one finds the universal ratio 
C to be 5.55878 -I- 0{1/N). To the leading order, the 
scaling function does not depend on the exact nature of 
the short-ranged interactions. Any order parameter that 
lifts the valley degeneracy from the zeroth LL will lead 
to identical scaling behavior. However, 1/A'' corrections 
may depend on the exact nature of the order parameters 

MM- 

Pseudo-magnetic fields may arise from specific de- 
formation of the graphene sheet, thus field profile is 
typically inhomogeneous. Hence it natural to study the 
behavior of the collection of the interacting fermions on 
honeycomb lattice when the fields are inhomogeneous. 
A recent numerical study [31] in the presence of a 
nonuniform, but real magnetic field, established that the 
catalysis mechanism survives even when the condition 
of uniformity of the field is relaxed. System finds itself 
in an ordered phase at sufficiently weak interaction by 
developing a local expectation value of the order param- 
eter, which more or less is found to follow the profile 
of the magnetic field. System finds itself in an ordered 
phase by lifting the chiral symmetry. A similar behavior 
of the local order parameter on the local strength of 
the field is been observed, when only a inhomogeneous 
pseudo-magnetic field penetrates through the system 
|3a |. The order parameter in presence of fictitious field 
breaks the TRS [s^. Hence, in presence of inhomo- 
geneous fields, we will expect the order parameter to 
develop a space modulated expectation value. However, 
we leave this issue for further investigation. 



V. SUMMARY AND DISCUSSION 

To summarize, we consider Dirac fermions in two 
spatial dimensions in the presence of real and pseudo- 
magnetic fields. For simplicity, we assumed both of 
them to be uniform and the real one is directed perpen- 
dicular to the graphene plane. Diagonalizing the free 
Hamiltonian under that circumstance we found that all 
the LLs at finite energy lack valley degeneracy, whereas 
the one at zero energy enjoys it. Taking into account 
the onsite and nearest-neighbor interaction we have 
shown that, the zero energy level splits by developing 
either a CDW or ferrimagnetic order when the real 
magnetic field is stronger than the fictitious one. Upon 
including the Zeeman splitting, the spin degeneracy 
from all the LLs is lifted and one expects to see plateaus 
in the Hall conductivity at all the integer factors of 
e'^/h. Previously the plateaus at odd integer fillings 
/ = ±3, ±5, • • • could not be observed due to the valley 
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degeneracy of the higher LLs. Here we proposed that 
a finite pseudo flux can, in principle, release the LLs 
from the valley degeneracy. However, due to distinct 
degeneracies of the LLs localized in the vicinity of two 
Dirac points, the plateaus in Hall conductivity at integer 
values of e^/ft., would be observed at non-integer fillings. 
A scaling behavior of the interaction induced gap at 
/ = 1 is presented, assuming either onsite Hubbard or 
nearest-neighbor repulsion is the dominant finite range 
interaction. Within the mean field approximation, the 
size of the gap is found to be enhanced in the presence 
of a finite pseudo flux. One way to test the enhancement 
of the gap is as follows. As mentioned in Ref. [T^ . 
the bulging of the graphene sheet does not take place 
everywhere on the sample. Only some portion of the 
flake develops a bump and experiences a flnite pseudo 
field. Therefore placing the system in a magnetic field 
one can measure the order parameter in two different 
domains of the flake, where the flake is bulged and the 
other where it is reasonably flat. This way one can 
determine the effect of the flnite pseudo flux on the 
interaction driven gap. 

In a finite system, pseudo field can be achieved 
by means of an inhomogeneous strain. Particular 
deformations that might give rise to such field was 



One feature 



proposed in Ref. [lOj and in Ref. 
of such deformation is that it leaves the lattice only 
with a C3 symmetry, which has also been confirmed 
experimentally As long as the extension of the 

bump in graphene flake and the magnetic length is 
larger than the lattice spaciiig, continuum description is 
valid. In recent experiment [l2|. It, sa 35 A at 6 = 350 T 
and the extension of the bump is '-^ 40 — 100^. Hence, 
both the length scales are much larger than the lattice 
spacing in graphene (a « 2.5A). However, for true 
thermodynamic limit, one requires the extension of the 
bump to be much larger than the magnetic length. In 
present experimental situation these two length scales 
are comparable. However, with a smoother bump with 
broader extension, one can get to the regime where the 
constrain is easily satisfied. 

Finally let us turn to the situation when b > B. 
Recently a uniform pseudo field of strength 350 T was 
produced by bulges in the graphene layer deposited on 
Platinum substrate [l2j . Hence this limit seems quite 
achievable. Particularly for B ^ the zero energy states 
near two Dirac points live on the same sublattice. Hence, 
an infinitesimal next-nearest-neighbor Coulomb repul- 
sion may open up a gap by spontaneously breaking time 
reversal symmetry (TRS) [36i|. The order parameter 
in the insulating phase is given by (4'J.(^2 <8) 17172)^(7) ■ 
The TRS is represented by an anti-unitary operator 
It = UtK, where Ut is unitary and K is complex 
conjugation operator. In the graphene representation 
Ut is given by 17175 = ci 65 /2 [lal- The correlated 
ground state under this circumstance corresponds to 



an intra sublattice Haldane circulating current, prop- 
agating in opposite direction on the two sublattices 
[37|. This state is named the quantum anomalous Hall 
insulator (QAH). The exact nature of the correlated 
ground state in presence of a pseudo-magnetic field 
is yet to be determined and we leave this issue for 
future investigation. However, in the absence of any 
gauge potentials {B — 0,b = 0) fluctuations preempt the 
appearance of QAH and stabilizes the spin Hall insulator 
(QSH) ground state. The correlated ground state in the 
QSH phase, breaks the TRS for each spin component 
and the order parameter reads as (5'J(t (8) «7i72)5'ct) [ssj . 

Here we focused only on the flnite-ranged components 
of the Coulomb interaction. On the other hand, its 
long range 1/r tail is unscreened due to the vanishing 
density of states at the charge neutral point. If one 
turns off the pseudo-magnetic fleld the long range 
Coulomb interaction will lead to a similar splitting 
of the LLs. However, it introduces an energy scale 
e^/elB, which is on the same order of the LL energy 
with commonly assumed dielectric constant e sa 5 
for Si02 substrate. This immediately contradicts the 
experimentally observed gap of ~ 100 K at / = 1, an 
order of magnitude smaller than the LL energy [28|. 
This issue can be resolved by assuming an order of 
magnitude larger value of the dielectric constant, which 
attributes the accumulation of water layer between the 
graphene flake and 5*102 substrate [1^. In recent work 
|40|, it was shown that even though the long range 
Coulomb interaction liberates all the LLs from the valley 
degeneracy, the critical strength of the disorder above 
which the / = 3 plateau is destroyed is much smaller 
than that for f — 1. This scenario, may explain the 
absence of odd integer Hall plateaus in graphene 0, HI] . 
In contrast, the presence of a flnite flux of pseudo field 
leads to an activation gap ~ 200 — 400 K for / = ±3, ±5 
Hall states (Table I). Consequently, the Hall states 
at those fillings should be much more robust against 
disorder. Thus, one may expect to see the quantization 
of Hall conductivity at odd integer factors of e^/h in a 
cleaner system. Assuming a sufficiently large dielectric 
constant, one can restrict oneself to the short-ranged 
components of the Coulomb interaction. Effect of weak 
long-ranged Coulomb interaction on the interaction 
mediated gap is discussed in Ref. [30| . 
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VII. APPENDIX 

Here we determine the quantization of the HaU con- 
ductivity that we mentioned at the end of the Sec. II 
and Sec. III. One notices that upon introducing a fi- 
nite pseudo-magnetic flux, aU the LLs do not enjoy equal 
degeneracies. LLs localized in the neighborhood of the 
Dirac point at K, have a degeneracy {B -f &)/27r per unit 
area, whereas those living in the vicinity of the other 
Dirac point, at —K, carry {B — 6)/27r states per unit 
area. Using the Stfeda formula, one can write the expes- 
sion of the Hall conductivity as 

after setting e — c — 1 [4l|. Here N is the electronic 
density in the bulk and the derivative with respect to 
B is taken at fixed chemical potential (^). The chemi- 
cal potential is measured from the charge neutral point. 
Moreover, all the LLs have additional two fold spin de- 
generacy. Upon changing the magnetic field -B, number 



of states below the chemical potential is changed to 

5N = flf+SB + nf_6B. (18) 

/-I- and /_ counts the number of LLs below the chem- 
ical potential (/i), but above the charge neutral point, 
with degeneracies {B ± b) /27r per unit area for each spin 
species, respectively. Therefore the Hall conductivity is 

^.y ^f+ + f-=f- (19) 

Hence, quantization of the Hall conductivity only takes 
the number of filled LLs below the chemical potential 
into account. That leads to the the announced result 
of the Hall conductivity in presence of real and pseudo- 
magnetic fields. However, one should notice that due 
to the different degeneracies of the LLs, quantization 
of the Hall conductivity at integer values of e'^/h is 
not associated with integer fillings. The filling factors 
explicitly depend on the ratio b/B. 
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